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ABSTRACT 

Context. Partially ionized plasma is usually described by single-fluid approach, where the ion-neutral collision effects are expressed 
by Cowling conductivity in the induction equation. However, the single-fluid approach is not valid for the time-scales less than ion- 
neutral collision time. For these time-scales the two-fluid description is better approximation. 

Aims. To derive the dynamics of magnetohydrodynamic waves in two-fluid partially ionized plasmas and to compare the results with 
those obtained under single-fluid description. 

Methods. Two-fluid magnetohydrodynamic equations are used, where ion-electron plasma and neutral particles are considered as 
separate fluids. Dispersion relations of linear magnetohydrodynamic waves are derived for simplest case of homogeneous medium. 
Frequencies and damping rates of waves are obtained for different parameters of background plasma. 

Results. We found that two- and single-fluid descriptions give similar results for low frequency waves. However, the dynamics of 
MHD waves in two-fluid approach is significantly changed when the wave frequency becomes comparable or higher than ion-neutral 
collision frequency. Alfven and fast magneto-acoustic waves attain their maximum damping rate at particular frequencies (for exam- 
ple, the peak frequency equals 2.5 ion-neutral collision frequency for 50 % of neutral Hydrogen) in wave spectrum. The damping 
rates are reduced for higher frequency waves. The new mode of slow magneto-acoustic wave appears for higher frequency branch, 
which is connected to neutral hydrogen fluid. 

Conclusions. The single-fluid approach perfectly deals with slow processes in partially ionized plasmas, but fails for time-scales 
smaller than ion-neutral collision time. Therefore, two-fluid approximation should be used for the description of relatively fast pro- 
cesses. Some results of single-fluid description, for example the damping of high-frequency Alfven waves in the solar chromosphere 
due to ion-neutral collisions, should be revised in future. 
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' 1. Introduction 

Astrophysical plasmas often are partially ionized. Neutral atoms 
, may change the plasma dynamics due to collisions with charged 
' particles. The ion-neutral collisions may lead to different new 
I phenomena in plasma, for example the damping of magnetohy- 
. drodynamic (MHD) waves (Khodachenko el al. 120041 Forteza 
■ et al. l2007l l. Solar photosphere, chromosphere and prominences 
[ contain significant amount of neutral atoms, therefore the com- 
. plete description of plasma processes requires the consideration 
' of partial ionization effects. 

Braginskii ( 1 19651 1 gave the basic principles of transport pro- 
cesses in plasma including the effects of partial ionization. 
Since this review, numerous papers addressed the problem of 
partial ionization in the different regions of solar atmosphere. 
Khodachenko and Zaitsev (I2002I I studied the formation of mag- 
netic flux tube in a converging flow of solar photosphere, while 
Vranjes et al. (2008) studied the Alfven waves in weakly ionized 
photospheric plasma. Leake and Arber ( I2005I I and Arber et al. 
(12007b studied the effect of partially ionized plasma on emerging 
magnetic flux tubes and concluded that the chromospheric neu- 
trals may transform the magnetic tube into force-free configura- 
tion. Haerendel ( |1992| l. De Pontieu and Haerendel ( 119981 ). James 
and Erdelyi ( I2002I I. James et al. ( I2004l i considered the damp- 



ing of Alfven waves due to ion-neutral collision as a mecha- 
nism of spicule formation. Khodachenko el al. (I2004I I and Leake 
et al. ( I2OO6I 1 studied the importance of ion-neutral collisions in 
damping of MHD waves in the chromosphere and prominences. 
Forteza et al. ( l2007l |2^, Soler et al. ( I2009al I2009bl l20T0b 
and Carbonell et al. (120 lOj studied the damping of MHD waves 
in partially ionized prominence plasma with and without plasma 
flow. 

All these papers considered the single-fluid MHD approach, 
when inertial terms in the momentum equation of relative ve- 
locity between ions and neutrals are neglected. The partially 
ionized plasma effects are described by generalized Ohm's law 
with Cowling conductivity, which leads to the modified induc- 
tion equation (Khodachenko el al. I2004I I. Ambipolar diffusion 
is more pronounced during the transverse motion of plasma 
with regards to magnetic field, therefore the Alfven and fast 
magneto-acoustic wave are more efficiently damped. The slow 
magneto-acoustic waves are weakly damped in the low plasma 
beta case. Moreover, Forteza et al. (I2007I I found that the damp- 
ing rate of slow magneto-acoustic waves derived through nor- 
mal mode analysis is different from that estimated by Braginskii 
( 119651 ). The problem of discrepancy between normal mode anal- 
ysis (Forteza et al. 120071 ) and energy consideration (Braginskii 
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11965b is still an open question and the present study attempts to 
shed light on it. 

The single-fluid approach has been shown to be valid for 
the time-scales which are larger than ion-neutral collision time. 
However, the approximation fails for the shorter time scales, 
therefore the two-fluid approximation, which means the treat- 
ment of ion-electron and neutral gases as separate fluids, should 
be considered. The two-fluid approximation is valid for the time- 
scales larger than ion-electron collision time, which is signifi- 
cantly short due to Coulomb collision between ions and elec- 
trons. 

In this paper, we study MHD waves in two-fluid partially 
ionized plasma. The particular attention is paid to the wave 
damping due to ion-neutral collisions and comparison between 
the wave dynamics in single and two-fluid approximations. We 
derive the two-fluid MHD equations from initial three-fluid 
equations and solve the linearized equations in the simplest case 
of a homogeneous plasma. 

2. Main equations 

We aim to study partially ionized plasma, which consists in 
electrons, ions and neutral atoms. We suppose that each sort of 
spaces has Maxwell velocity distribution, therefore they can be 
described as separate fluids. Below we first write the equations 
in three-fluid description and then perform consequent transition 
to two-fluid and single-fluid approaches. 



2.1. Three-fluid equations 

The fluid equations for each spaces can be derived from 
Boltzmann kinetic equations and they have the following forms 
(Braginskii [T965l Goedbloed & PoedtsEOQl 



dt 



+ V ■ (n,y,) = 0, 



5 + v.(«,y,) = o, 

at 

^ + v.(«,.y„) = o, 

at 



I dV- 

_i+(y..V)y,. 



(1) 
(2) 
(3) 



where nia, ria, Pa, Ta, Va are the mass, the density, the pressure, 
the temperature and the velocity of particles a, E is the electric 
field, B is the magnetic field strength, qg is the heat flux density 
of particles a, /?„ is the change of impulse of particles a due to 
collisions with other sort of particles, Qa is the heat production 
due to collisions of particles a with other sort of particles, 7r„ is 
the off-diagonal pressure tensor of particles a, e = 4.8 x 10""^ 
statcoul is the electron charge, c = 2.9979 x 10'° cm s"' is the 
speed of light and k - 1.38 x 10 erg K ' is the Boltzmann 
constant. The double dot indicates that a double sum over the 
Cartesian components is to be taken. Plasma is supposed to be 
quasi-neutral, which means rie - Zn\. In what follows we con- 
sider hydrogen ions and hydrogen neutral atoms which imply 
Z - \. The description of the system is completed by Maxwell 
equations which have the forms (without displacement current) 



VxE^ , 

c dt 

„ 47r . 
V X B = — 7, 

c 

where 

j = -eiieiVe - y,) = -e«f. 



(11) 



(12) 



(13) 



is the current density. 

In the case of Maxwell distribution in each sort of particles, 
Ra and Qa are expressed as the following (Braginskii 1 1 965l l : 



Re 
Ri 

R„ 

Qe 

Qi- 

Qn 



■■ -a„(n - y,) - a,„(Ve - y„), 
-ff,v(y,- - y,) - a,„(y,- - y„), 

: -aneiVn - V e) ' ani(V„ - Vi), 
■ a.iiV, - ViWe + a,„{Ve - V„)Ve, 
aiAVi - V,)Vi + ai„{Vi - V„)Vi, 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



aha are coefficients of friction between particles a 



where Oab 
and b. 

For time scales longer than ion-electron collision time, the 
electron and ion gases can be considered as a single fluid. This 
_ _ 1 1 r.\ , significantly simplifies the equations taking into account the 

' ' I c ' I f V smallness of electron mass with regards to the masses of ion and 

neutral atoms. Then the three-fluid description can be changed 
]^ \ by two-fluid description, where one component is ion-electron 

-V pi-V -Tii+Zeni \E^ — y, x B \+Ri, (5)gas and the second component is the gas of neutral atoms. 



IdV 

m„n„ + iV„ ■ V)y„ I = -Vp„ -V ■7T„ + R„ 



(6) 



2.2. Two-fiuid equations 



Summing of Eqs. © and (|5]), Eqs. O and ([8]) and first two equa- 
. tions of Eq. ( fTOb . we obtain (after neglecting the electron inertia 

fi^ki^^ + (Ve ■ V)r. I -i-p.V -Ve + TTe ' , = -V ■ Op + Qe(1) ™d the viscosity cff'ect expressed by off'-diagonal pressure tensor 



3 IBT \ 

^riik {-g^+(Vi- V)r,J + p,V ■ Vi + 71 i : Vy,- = -V ■ qi + Q, (8) ^7 

^«„^j^ + (y„ ■ v)r„J+p„v-y„+7r„ : vy„ = -v-^„ + e„(9) 



+ V ■ (n,y,) = 0, 



+ V ■ («„y„) = 0, 



(20) 
(21) 



Pe - iiekTe, pi - iiikTi, p„ - n„kT„, 



(10) 



mmi I ^ + (Vi ■ V)y,) = -Vpi, + -jxB + ^j- 
ot / c erie 



(a,„ + ae„)(Vi 
IdV, 



dt 



+ {V„ ■ V)V„ 



dt 



+(VrV)Pie+7Pie'^-Vi 
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(22) 3. Linear MHD waves 

,„ We consider the simplest case of static and homogeneous plasma 

-J+(«in+Q'«!)(^^/-V'„), (23) ^jj-jj homogeneous unperturbed magnetic field. Then the lin- 
earized two-fluid equations followed from Eqs. (l20t-(l25Tl and 
Oei .7 , _ .. .. (l27Ti are (neglecting the Hall term and the collision between neu- 



(r-i)^/+(r-iM«(V,-y„)-v,+ 



(y-l)a,„(V,-V„yV,+ ^ ' ^^P' +yp,V-^-(y-l)V-(g,+g,). (24) ^ +p/oV ■ v 



trals and electrons i.e. Ue,, <*c a,-,,): 
0, 



dpn 
dt 



+ {V„ ■ V)p„ + yp„V ■ V„ = -(r - l)ai„(V/ 



(y - l)a,„(V„ - V,) ■ V„ - (r - 1)V ■ q„. 



V„) ■ V„+ 
(25) 



where pie - p, + pe is the pressure of ion-electron gas and y - 
CplCv = 5/3 is the ratio of specific heats. 

The Ohm's law is obtained from the electron equation (Eq.|4|i 
after neglecting the electron inertia (i.e. the left-hand side terms) 
and it has the form 



E+-ViXB+—Vpe 

c en^ 



— Tl—J (V/-K)+ JXB.(26) 



Maxwell equation (Eq. fTTI) and Ohm's law (Eq. 
the induction equation 



dt 



Vx 



ca,,„(Vi - V„) 



cVpe 



■ Vx (t/V X B) - Vx 



lead to 



jxB 



where 



'7 = 



C-(aei + Sen) 



Ana- 



(27) 



(28) 



is the coefficient of magnetic diffusion. 

The coefficient of friction between ions and neutrals (in the 
case of same temperature) is calculated as (Braginskii ll965b 



n,n„OT,„cr,„: 



(29) 



where ;7z,„ = OT;m„/(m/ + m„) is reduced mass and cr,„ - ji{ri + 
r„)^ - Anrj is the collision cross section between ions and neu- 
trals. 

Collision frequency between ions and neutrals is then 



nin„mi, 



niirii + lllnVln 

, In [kT 

3 «; + «« V mi 



3 m,«/ -I- m„n„ 



32nrf 



5 X 10 



-12 



n, + n„ 



(30) 



where the atomic cross section nrf = 8.7974 x 10"'^ cm^ is used 
and T is normalised by 1 K. The chromospheric temperature of 
10"* K and hydrogen ion and neutral number densities of 2.3 x 
10'° cm-^ and 1.2x10'° cm^^ (Fontenla et al. [19901 model FAL- 
3) give the collision frequency as 4 s"'. 

For time scales longer than ion-neutral colUsion time (l/v,„), 
the system can be considered as a single fluid (the full equations 
of single-fluid MHD including neutral hydrogen are presented in 
AppendixA). However, when the time scales are near or shorter 
than ion-neutral collision time, then the single-fluid description 
is not valid and the two-fluid equations should be considered. 
In what follows we study the linear MHD waves in two-fluid 
description. 



dt 
dt 



+ P«oV ■ v„ = 0, 



(31) 
(32) 



P,o5 = -"^Pie - l^V(Bo ■ b) + ^(Bo ■ V)Z- + X b- 

dt An An Anerie 



CtiniVi - v„), 

P 0— = -^p' - 

dt " Anerif 



V X ft + Qf,„(V; - V„), 



^ = (Bo ■ V)v,- - BoV ■ V,- + T]V^b + ^Vx(v,- 
dt eUe 



v„), 



V, = 0, 



0, 



(33) 
(34) 

(35) 

(36) 

(37) 



where (p',) are perturbations of ion (neutral) density, v, (v„) 
are the perturbations of ion (neutral) velocity, p'.^ {p'^) are the 
perturbations of ion-electron (neutral) gas pressure, b is the 
perturbation of the magnetic field, and p,o - miniQ,p„Q - 
m„n„Q, Pie, p„, Bq are their unperturbed values, respectively. Eqs. 
(I3ni-(l32li and Eqs. (I36]l-(l37li lead to the expressions 



/ 111 1 I 

Pie = C„p,., p„ = C,„p„, 



(38) 



where c^,- = yfrpielpKi and c,,, = yfrpiTpi^ are sound speeds of 
ion-electron and neutral gases, respectively. 

Below we consider the unperturbed magnetic field, B , di- 
rected along the z axis and wave propagation in xz plane i.e. 
djdy - 0. Then Eqs. (I3ni-(l37b can be split into Alfven and 
magneto-acoustic waves. 

3. 7 . Alfven waves 

Let us assume the Alfven waves polarized along y axis. We in- 
tend to study the damping of Alfven waves due to collision be- 
tween ions and neutrals. Therefore, we neglect the magnetic dif- 
fusion for simplicity. Then, Eqs. OTli-dJTli give 



dViy 

'df 



B- 



db. 



Anpio dz Pit) 



(viy - v„v), 



dv„r Oin ^ ^ 

dt p„o ■ 



db, 
'dt 



B 



dViy 

~dl' 



(39) 
(40) 
(41) 



Fourier analyses assuming disturbances to be proportional to 
exp{i{k.z - wf)] give the dispersion relation 



a£,i£,nTn + izn — a^nVJ — i — 0, 



(42) 
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^ =0.5 ^ =0.5, k =0 




5 10 15 2 4 6 8 10 




Fig. 1. Real (nr^) and imaginary {mi) parts of Alfven wave fre- 
quency, vj, vs a (where a - k-VA/vi„). The blue line corresponds 
to the solution of single-fluid dispersion relation, i.e. Eq. (l44l i 
and red asterisks are the solutions of two-fluid dispersion rela- 
tion, Eq. (l42l i. The values are calculated for 50% of neutral hy- 
drogen, ^,,=0.5. 
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Fig.2. Same as on Fig.l, but for 10% of neutral hydrogen, 
f„=0.1. 



where 



a> 

m = , a 

k-VA 



Oh 



Kva . Pi . Pn B, 

, = — , = — , Va = . 

'^in Pi) po ^j4-npo 



Vin - , Po = P/0 +P;iO- 

PO 



(43) 



The same dispersion relation can be obtained from linear 
single-fluid equations retaining the inertial term in Eq. ( IA.6I 1. 
The dispersion relation of the Alfven waves in linear single-fluid 
equations without the inertial term can be easily derived as 



m" + ia^„m —1=0. 



(44) 



Fig. 3. Frequency and damping rate of diff'erent wave modes 
in two-fluid MHD vs a (where a - kvA/vj,,). Frequencies and 
damping rates are normalized by kvA- Red asterisks correspond 
to fast magneto-acoustic mode and green diamonds correspond 
to usual slow magneto-acoustic mode. The mode with the blue 
squares is the new sort of slow magneto-acoustic wave ("neutral" 
slow mode), which arises for larger wave numbers. This mode 
has only imaginary frequency for small wave numbers, which is 
not shown in this figure. Frequencies and damping rates are cal- 
culated for the waves propagating along the magnetic field. The 
neutral hydrogen is taken to be 50% (^„=0.5). Here we consider 

Csn/VA = 0.5. 



The solution of Eq. (l44l l is 



which for a^l < 2 gives the damping rate 



2a)i 



(45) 



(46) 



in full coincidence with Braginskii ( 119651 ). On the other hand, 
the condition a^j^ > 2 in Eq. ( |46] | retains only imaginary part, 
which gives the cut-off wave number 



kr^ 



2Vm 



(47) 



The value of cut-off' wave number has been obtained recently 
by Barcelo et al. ( 120101 ). Hence, the waves with higher wave 
number than k^ are evanescent. However, it might be incomplete 
conclusion as the complete treatment requires inclusion of in- 
ertial terms, and therefore dealing with Eq. (l42l) instead of Eq. 
( l44l i. The first term in Eq. (l42l i is important for the high fre- 
quency part of wave spectrum and could not be neglected. We 
demonstrate it by solutions of Eqs. (l44l ) and ( l42b . 

Fig. 1 displays the solutions of single-fluid (Eq. |44l blue 
lines) and two-fluid (Eq. |42l red asterisks) dispersion relations 
for ^„ - 0.5. We see that the frequencies and damping rates of 
Alfven waves are same in single-fluid and two-fluid approaches 
for low-frequency branch of spectrum (small a). But the behav- 
ior is dramatically changed when the wave frequency becomes 
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comparable or higher than the ion-neutral collision frequency, 
v,„, i.e. for a > 1 . The damping time linearly increases with a and 
the wave frequency becomes zero at some point in single-fluid 
case (blue lines). The point where the wave frequency becomes 
zero corresponds to the cut-off wave number kc of Barcelo et al. 
(12010b . However, there is no cut-off wave number in solutions 
of two-fluid dispersion relation (red asterisks): Eq. (l42l) always 
has a solution with a real part. Therefore, the occurrence of the 
cut-off wave number in single-fluid description is the result of 
neglecting the inertial terms in the momentum equation of rela- 
tive velocity between ions and neutrals. Therefore, Eq. ( 142b is the 
correct dispersion relation for the whole spectrum of waves. But 
the dispersion relation (l44l) is still good approximation for lower 
frequency part of spectrum. Another interesting point of two- 
fluid approach is that the damping rate (i.e. cjj) attains its max- 
imal value at some wave-lengths for which a x 2.5. The damp- 
ing rate decreases for smaller and larger a. This means that the 
waves, which have the frequency in the interval v,„ < w < 10 v,„, 
have stronger damping than other harmonics of spectrum. This 
is totally different from single-fluid solutions, which show the 
linear increase of damping rate with increasing wave number 
(lower panel, blue line). 

Fig. 2 displays the same solutions as on the Fig. 1, but for 
^„ - 0.1. The solutions have basically same properties as those 
with ^„ = 0.5. However, the wave length with maximal damping 
rate is now shifted to a ^ 10. 



3.2. Magneto-acoustic waves 

Now let us turn to magneto-acoustic waves. We consider the 
waves and wave vectors polarized in xz plane. Then Eqs. (ISTb - 
dJTl i are written as (magnetic diffusion is again neglected) 



dp- (dvi, dvA 



ldv„ 



dVn: 



dt ^"^ \ dx dz 



dv 



1 dp'ie 



= 0, 

B. db. 



(48) 
(49) 



B- db^ a in 
dt pio ox 4-npiQ dx 4-npio dz p,o 



1 dpi am 

^ (Viz - v„J, 

dt pio dz Pio 



dvi- 
~dt 

dv. 



1 dp„ ai„ 
dt pno dx p„o 



dt 

db, 
'dT 



j_d_p[^ 

PnO dz 

dVir 



PnO 



-{Viz - V'nz), 



dt 
dt 



■ dz 



■ ypu 



+ ypn 



I dvi_^ 


dvi. 


\ dx 




'dv„x 


dVnz 


dx 


"IF 



= 0. 



(51) 
(52) 
(53) 
(54) 
(55) 
(56) 



Fourier analysis with exp{i{kxX-\-kzZ— dit)^ and some algebra 
give the dispersion relation 

vloj [cJ - k\c% + c]„^„ + Vl)a? + (cl^i + cl,^„)kXvl] - 




Fig. 4. Damping rate of fast (upper panel) and slow (lower panel) 
magneto-acoustic waves, i.e imaginary part of cj normalized by 
kvA, vs a (where a - A:v/i/v,„). The blue solid lines corresponds 
to the solution of single-fluid dispersion relation and dashed 
line corresponds to the slow magneto-acoustic damping rate of 
Braginskii (the expressions are used from Forteza et al. 12007b . 
Red asterisks are the solutions of two-fluid dispersion relation, 
Eq. (|57] |. The values are calculated for 10% of neutral hydrogen, 
^,,=0.1, and for Cs„/va =0.1. The damping rates are calculated 
for the waves propagate with 45" angle with regards to the mag- 
netic field. 



^i^lojico^ - cik^) yao' - k^vW + clk\klvl - ^iU?-)\ - 



■2, ,2 



•(^«-2)w" 



■clk^oj^ 



■klvl))]^0, (57) 



where k = ^k^ H- k^. 

The dispersion relation (ISTT i is seventh order equation with 
o), therefore it has 7 different solutions. For smaller wave num- 
bers (or smaller frequencies) four of the solutions represent the 
usual magneto-acoustic waves, while 3 other solutions are purely 
imaginary and are probably connected to the vortex modes (with 
Re((D) = 0) damped due to ion neutral collisions. The vortex 
modes are solutions of fluid equations and they correspond to the 
fluid vorticity. The vortex modes have zero frequency in the ideal 
fluid, but may gain purely imaginary frequency if dissipative 
processes are evolved. The two vortex modes are transformed 
into oscillatory modes for shorter wavelengths (see the next 
paragraph). Then we have two fast magneto-acoustic modes, 
four slow magneto- acoustic modes and one vortex solution with 
purely imaginary part. In what follows, we consider that the 
temperatures of all three species are equal i.e. T, - Tg - T„, 
which gives c^. = jpielpio = y(pi + Pe)lpm = yk(Ti + Te)lmi = 
2ykT„/m„ = 2yp„/p„o = 2c^„. 

Fig. 3 displays all oscillatory solutions of two-fluid disper- 
sion relation for ^„ - 0.5 (only the modes with positive frequen- 
cies are shown). The wave propagation is parallel to the mag- 
netic field and we use Cs„/va = 0.5, where Cj,, is the sound speed 
of neutral hydrogen. For smaller wave numbers, A; < 3.5 Vi„/vA, 
there are two usual magneto-acoustic modes, fast (red asterisks) 
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Fig. 5. The same as in Fig. 4 but for the parallel propagation, i.e. 
k, = 0. 



and slow (green diamonds). However, for larger wave-numbers, 
k > 3.5 Vi„/vA one additional sort of slow magneto-acoustic 
mode with strong damping rate (blue squares) arises. The "neu- 
tral" slow mode is connected with neutral atoms. For higher 
frequency range, i.e. for those larger than ion-neutral collision 
frequency, neutral gas does not feel the ions, therefore it sup- 
ports the propagation of additional oscillatory wave mode. This 
mode obviously disappears for lower frequency as the collisions 
couple ions and neutrals and they behave as a single fluid. In 
other words, for lower frequencies (or small wave numbers) this 
mode has zero real part, but non-zero imaginary part (not shown 
in the figure). Therefore, more correct statement is that the os- 
cillatory mode transforms into non-oscillatory vortex mode for 
smaller wave numbers. The fast magneto-acoustic modes decou- 
ple from the slow waves for the parallel propagation and they 
have the same behavior as Alfven waves. Therefore, the plot of 
fast magneto-acoustic waves is similar to that of Alfven waves 
(see Fig.l). 

It is useful to compare the solutions of two-fluid disper- 
sion relation to those obtained in the single-fluid approach. 
The damping of fast and slow magneto-acoustic waves have 
been derived from the energy equation by Braginskii (|1965l l. 
Khodachenko et al. (|2004] |, Khodachenko & Rucker (2005) 
and through normal mode analysis by Forteza et al. (2007|. 
Damping rates are same in both considerations for fast magneto- 
acoustic waves, but they disagree for slow magneto-acoustic 
waves (Forteza et al. I2007I I. Namely, slow magneto-acoustic 
waves have damping for purely parallel propagation in the case 
of Braginskii, while the damping is absent in the case of Forteza 
et al. Fig. 4 and 5 show the damping rates of fast and slow 
magneto-acoustic waves vs a (i.e. k) for the propagation angles 
of 45'^ and O", respectively. Red asterisks are the solutions of 
two-fluid dispersion relation- Eq. ( |57] |. The blue solid lines cor- 
respond to the solutions of single-fluid dispersion relation from 
Forteza et al. ( I2007l i. The dashed line corresponds to the slow 
magneto-acoustic damping rate of Braginskii ( |1965t . Here we 
use Cs„/\'A = 0. 1 so the plasma /3 is small enough. Neutral hydro- 
gen concentration is taken to be 10%. The fast magneto-acoustic 
waves have essentially same dynamics as the Alfven waves. For 
the smaller wave numbers (or smaller frequencies) the two-fluid 



Fig. 6. Damping rate of fast (upper panel) and slow (lower panel) 
magneto-acoustic waves vs a . Red asterisks are the solutions of 
two-fluid dispersion relation. Blue lines correspond to the solu- 
tions of Braginskii. The values are calculated for 50% of neutral 
hydrogen and Cs„/va =0.1. The damping rates are calculated for 
the waves propagate along the magnetic field. 



and single-fluid approaches give the same results, but for the 
larger wave numbers the damping rate is decreased in two-fluid 
description as in the case of Alfven waves. On the other hand, 
the slow magneto-acoustic waves have similar damping rates in 
both approaches. There is a small discrepancy between damping 
rates for the waves propagating with 45° degree about the mag- 
netic field, but all the three cases (two-fluid waves, single-fluid 
waves and energy consideration) give similar results. The paral- 
lel propagation reveals an interesting result: two-fluid solutions 
are exactly the same as those obtained by Braginskii (the damp- 
ing rate obtained by Forteza et al. 120071 is zero for the parallel 
propagation). Therefore, the discrepancy between damping rates 
of slow magneto-acoustic waves obtained by Forteza et al. 120071) 
and Braginskii (119651 ) is again caused by neglecting the inertial 
terms in the momentum equation of relative velocity (Eq. lA.6b . 
Braginskii ( 1 19651 ) used energy equation to calculate the damping 
rate, therefore his solution agrees to that obtained in two-fluid 
approach. 

Fig. 6 shows the comparison of damping rates in two-fluid 
approach and those obtained by Braginskii ( 119651 ) in the case 
of ^„ - 0.5, Csn/vA = 0.1 and parallel propagation. The fast 
magneto-acoustic waves have the same bihaviour as the Alfven 
waves (see lower panel of Fig. 1 ), which is significantly different 
from Braginskii (11965b and Forteza et al. (12007b . But, the slow 
magneto-acoustic waves have the same damping rate as that of 
Braginskii (119651 ). On the other hand, the damping rate of slow 
magneto-acoustic waves becomes different from the solution of 
Braginskii (11965b for higher plasma Fig. 7 shows the same 
as in the Fig. 6 but for c„,/v'a - 0.5. The damping rate of slow 
magneto-acoustic waves now begins to deviate from the solu- 
tion of Braginskii for higher wave numbers. The behavior of fast 
magneto-acoustic waves remains the same. 
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^ =0.5, k =0, fast waves 




Fig. 7. The same as in Fig. 6 but for Csnh'A - 0.5. 



4. Discussion 

Some parts of the solar atmosphere contain large number of 
neutral atoms: most of atoms are neutral at the photospheric 
level, but the ionization degree rapidly increases with height due 
to increased temperature. Solar prominences also contain neu- 
tral atoms. Neutral atoms may change the dynamics of plasma 
through collision with charged particles. For time scales longer 
than ion-neutral collision time, the partially ionized plasma can 
be considered as one fluid, because collisions between neutrals 
and charged particles lead to the rapid coupling of the two fluids. 
Then the equation of motion is written for the center of mass ve- 
locity, and the motion of species is considered as diffusion with a 
small velocity as compared to the velocity of center of mass. The 
corresponding collision terms appear in the equation of motion 
for the relative velocity (between ions and neutrals) and in the 
generalized Ohm's law. Neglecting the inertial term in the equa- 
tion of motion for the relative velocity one simplifies the equa- 
tions and traditional induction equation with Cowling conduc- 
tivity is obtained (Braginskii 119651 Khodachenko et al. I2004I I. 
The inertial terms (left-hand side terms in Eq. IA.6I 1 is smaller 
than the collision term (the last term in the same equation), but 
becomes comparable for time scales near ion-neutral collision 
time. Therefore, it can be neglected only for longer time-scales. 

On the other hand, for the time scales less than ion-neutral 
collision time, the both fluids may behave independently and the 
single-fluid approximation is not valid any more. Then the two- 
fluid approximation, when ion-electron and neutral atom gases 
are treated as separate fluids, should be considered when one 
tries to model the processes in partially ionized plasmas. 

The normal mode analysis of two-fluid partially ionized 
plasma shows that frequencies and damping rates of low- 
frequency MHD waves are in good coincidence with those found 
in the single-fluid approach. However, the waves with higher fre- 
quency than ion-neutral collision frequency have significantly 
different behavior. Alfven and fast magneto-acoustic waves have 
maximal damping rates in particular frequency interval peaking, 
for example, at the frequency u - 2.5 v,„, (v„, is the ion-neutral 
collision frequency) for ^„ = 0.5 and at the frequency o) - 10 v,„ 
for ^„ =0.1. The damping rates are reduced for higher frequency 
part of wave spectrum (note that the damping rates are linearly 



increased in the single-fluid approach). Therefore, the statements 
concerning the damping of high-frequency Alfven waves in the 
solar chromosphere due to ion-neutral collisions should be re- 
vised. Careful analysis is needed to study the damping of high 
frequency Alfven waves for realistic height profile of ionization 
degree in the chromosphere. 

Another important point concerning the Alfven waves in par- 
tially ionized plasma is the cut-off wave-number, which appears 
in the single-fluid approach (Barcelo et al. 120101 ). Barcelo et al. 
( 120101 ) found that the Alfven waves with larger wave numbers 
than the cut-off value are evanescent in partially ionized and re- 
sistive plasmas. However, our two-fluid analysis shows that there 
is no cut-off wave number due to ambipolar diffusion (see Fig. 1, 
red asterisks). Therefore, the appearance of cut-off wave number 
in the single-fluid approach is the result of neglecting of inertial 
term in the equation of motion for relative velocity. It is possible 
that the cut-off wave number due to usual magnetic resistivity 
is caused by neglecting of electron inertia, therefore the cut-off 
may completely disappear in three-fluid approach. The cut-off 
wave numbers also appear for fast magneto-acoustic waves in 
partially ionized and resistive plasmas (Barcelo et al. l2010l l. We 
suggest that this may also caused by the neglecting of inertial 
term. However, this point needs further study. 

The two-fluid approach reveals two different slow magneto- 
acoustic modes when the slow wave time scale becomes shorter 
than ion-neutral collision time (Fig. 3). The different slow modes 
correspond to ion-electron and neutral fluids. But, only one slow 
magneto-acoustic mode remains at lower frequency range as in 
usual single-fluid approach. This is easy to understand physi- 
cally. When the wave frequency is lower than ion-neutral col- 
lision frequency, then the two fluids are coupled due to colli- 
sions and only one slow magneto-acoustic wave appears. The 
mode connected with the neutral fluid ("neutral" slow mode) 
has only imaginary frequency in this range of wave spectrum 
(not shown in the figure). This means that any slow wave type 
changes (density, pressure) in the neutral fluid is damped faster 
than wave period due to collisions with ions. The "neutral" slow 
magneto-acoustic wave has similar properties as ion magneto- 
acoustic waves. 

It must be mentioned that the two-fluid approach of partially 
ionized plasma clarifies the uncertainty concerning the damping 
rate of slow magneto-acoustic waves found in the single-fluid 
approach. It was found that the normal mode analysis and en- 
ergy consideration method (used by Braginskii il965) lead to dif- 
ferent expressions for the damping rate of slow waves (Forteza 
et al. I2007I I. We found that the damping rate obtained in two- 
fluid approach is in good coincidence with the damping rate of 
Braginskii derived from the energy treatment (see lower panels 
of Fig. 5-6). Therefore, it seems that the discrepancy is again 
caused by neglecting the inertial terms in the equation of motion 
for the relative velocity in the single-fluid approach. Braginskii 
( 119651 ) used the general energy method for the estimation of 
damping rates, and probably this is the reason why his results 
agree to those found in two-fluid approach. 

Here we have considered only neutral hydrogen as a compo- 
nent of partially ionized plasma. However, other neutral atoms, 
for example neutral helium, may have important effects in MHD 
wave damping processes. Soler et al. ( 1201 Oi l made the first at- 
tempt to include the neutral helium in the single-fluid description 
of prominence plasma. They concluded that the neutral helium 
has not significant influence on the damping of MHD waves. 
However, two-fluid approach may give some more details about 
the effects of neutral Helium on MHD waves, therefore it is im- 
portant to study this point in the future. 
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5. Conclusions 

1. Frequencies and damping rates of low frequency MHD 
waves in the two-fluid description are similar to those 
obtained in the single-fluid approach. But high-frequency 
waves (with higher frequency than the ion-neutral collision 
frequency) have completely different behavior 

2. Alfven and fast magneto-acoustic waves have maximal 
damping rates at some frequency interval peaking at particu- 
lar frequency. The peak frequency is 2.5 v,„, where v,„ is the 
ion-neutral collision frequency, for 50% of neutral hydrogen. 
For 10% of neutral hydrogen, the peak frequency is shifted 
to 10 v,„. The damping rate is reduced for higher frequencies, 
therefore the damping of high-frequency Alfven waves in the 
solar chromosphere with realistic height profile of ionization 
degree needs to be revised in future. 

3. There are two types of slow magneto-acoustic waves in the 
high-frequency part of wave spectrum: one connected with 
the ion-electron fluid and another with the fluid of neutrals. 

4. There is no cut-off frequency of Alfven waves due to am- 
bipolar diffusion. The cut-off frequency found in the single- 
fluid approach is caused by neglecting the inertial terms in 
the momentum equation of relative velocity. 

5. The damping rate of slow magneto-acoustic waves is sim- 
ilar to Braginksii ( 119651 ) in low plasma approximation. 
The deviation from the Braginskii formula found by normal 
mode analysis in single-fluid approach (Forteza et al. I2007I I 
is probably caused by neglecting the inertial terms. 
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Appendix A: Single-fluid equations 

We use the total velocity (i.e. velocity of center of mass) 

PiVi +p„Vn 



Pi + Pn 

relative velocity 



(A.1) 



(A.2) 



and total density 

p=Pi+p„. (A.3) 
Eqs. (I20li-(|25]| and (|27] | lead to the system; 



^ + v-(py) = o, 

at 



(A.4) 



+ Piy ■ V)y = - + ij X B - V • (^i^npww), (A.5) 
ot c 

^ +{V ■V)w + (w ■ V)V + ^„{w ■ V)w - {w ■ V)^iW = 
ot 

+ — rJ X B + —J - ——w, (A.6) 



P^i P^n I Cp^i 

%+{V-V)p + ypV ■ V - ■ V)p - ypV ■ (f,H')+ 
at 

(w ■ V)pie + yp,>V -w ^{y- 1) 5-5—7" + (r - ^)anW - 

{y-l)—jW + — U-'^)Pe+yPe'^-^-{y-lW-{qi+qe+qnUf^-l) 

^=Vx(yxB) + Vx(£^)-Vx(,;VxB)-Vx(^] + 



Vx(^^Uvx(^„H'xB), 

\ erie j 



(A.8) 



where p = pe + pi + p„, ^/ = p,7p, ^„ = p„/p and a„ = a,-,, + Oe. 
The Ohm's law is now 

E+-VxB+ — Wpe - — T-^ — J w+ jxB-—wxB. {A.9) 

nLy^L 



c en,, 



e^nj eUf eerie c 



Neglecting the inertia terms i.e. all left hand side terms in 
Eq. ( IA.6I ) we have 

G . „ aen . 

w 1 J X B + J. 

a„ ca„ eriean 

Then the induction equation takes the form 

^ = Vx(y X B) + ^Vx ( YPiZI9\ _ vx {r^rV x B) - 

ot ^ \ / 

— Vxf ^~^^^" (Vxg)xg)-Vx(^Gxg) + 



Vx{-^{{V X B) X B) X B 

\Ana„ 



where e - ae„/a„, G - ^n^Pei - ^i^Pn and 



4-no- Ane-^nt a„ 



(A. 10) 



(A.ll) 



These equations are traditionally used for the description of 
partially ionized plasmas. 



